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ANALYTICAL CALCULATION OF (|Ci| 2 ) 


Consider two random gaussien vectors Ei = J2n=i c n 
and E’i = X^=i 4) where c n and c' n are random com¬ 
plex gaussien variable with (c n ) = (c' n ) = 0 and (|c n | 2 ) = 
(K| 2 ) = 1. Let us consider the correlation C\\ 
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Let us calculate (|Ci| 2 ) with the Law of large Numbers. 
We get: 
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Here, the n ^ p terms do not contribute to (|Ci| 2 ), be¬ 
cause c n , c^, c p and c' p are statistically independent. We 
get thus: 
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Since |c n | 2 and \c'n\ 2 are also statistically independent, 
we get: 
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This proves rigorously the result obtained by Monte 
Carlo for (|Ci| 2 ) . 


EXPERIMENTAL SETUP AND DATA ANALYSIS 
DETAILS 


Figure [T| shows the setup we used to study the open 
channels by measuring £ 2 ( 2 ;, //, t,p) and by calculating 
C 2 by : 
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FIG. 1: Experimental setup. E\\ incoming field; E 2 : out¬ 
going field; R( 1), R( 2): reference fields of polarisation p— 1 
and p— 2; PBS1, PBS2: polarized beam splitters; A/2: half 
wave plate to control £ 1 , R( 1) and R( 2) respective power; 
BS: beam splitter; M: mirror; DL: diffusing liquid; MOI and 
M02: microscope objectives; Q, Q’: camera plane, and cam¬ 
era conjugate plane with respect to M02; S: sample outgoing 
plane. P’ and P: MOI and M02 pupil planes. 


The setup consists of a Mach-Zehnder off-axis interfer¬ 
ometer with two orthogonally polarized reference beams. 
The light emitted by a A = 532 nm, 70 mW laser is 
split into a reference and an object field using a po¬ 
larizing beam splitter (PBS1). The studied sample is 
a ZnO powder slab with thickness l = 22 pm =b 7 pm 
deposited on a microscope cover slide. In order to maxi¬ 
mize the collection of both input and output modes, the 
sample is positioned between two microscope objectives: 
MOI (NA = 0.9 air, x60) in the powder side, and M02 
(NA = 1.4 oil, x60) in the cover slide side. 

A tank (1.5 cm thick) filled with viscous diffusing liq¬ 
uid (glycerol + concentrated milk) is positioned in front 
of MOI to randomize the illumination structure in both 
time and space. The incoming field E\ is therefore ran¬ 
domly distributed over all the incoming modes and varies 
in time. Thus, considering \t — t'\ > 100 ms, the fields 
£l(£) and £i(t') are uncorrelated. 

Measurement of the outgoing field £2 is holograph¬ 
ically performed. The two orthogonally polarized refer¬ 
ence beams R(p) (where p = 1, 2 is polarization) interfere 
with the outgoing fields E 2 (p ), and the interference pat¬ 
tern / = J2 P = 1 2 I R(p) + £ 2 (p)| 2 i s recorded on the CCD 
sensor (10 Hz, 1340 x 1040 pixels with Ax = 6.45 pm 
pitch). This configuration makes it possible to calculate 
from I the complex amplitudes E 2 (p) of the outgoing 
fields along both polarizations p = 1,2 directions by fil- 
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FIG. 2: |i7| 2 (a,d); x2 zoom of \H \ 2 near zone 1 (b); \E 2 R*\ 2 
calculated by Eq. [6] for p = 1 (c). In (c), the dashed squares 
of size M =724, 512, 362 and 256 pixels illustrate the ability 
to select different zones of the sample to calculate (|O 2 1 2 )- 
Images are obtained with ZnO sample (a-c), and without (d). 


tering, in the Fourier space, the desired +1 grating order 
i.e. E 2 {p)R*{p) Q. 

The Fourier spacial filtering is illustrated by Fig. [2j 
The 1340 x 1040 holograms were cropped to 1024 x 1024 
(the imaged area is therefore L = 77 fi m) and two 
frames holograms H were calculated from successive 
frames: i/(t 2n ) = /(^ 2 n) — /(£ 271 + 1 )• The Fourier holo¬ 
gram H(k x , k y ) = FFT [H (x, y)\ (where FFT is the Fast 
Fourier transform) is then calculated. \H \ 2 is displayed 
on Fig. [2fa). Four bright circular regions can be ob¬ 
served. They correspond to the reconstructed image of 
the M02 pupil. Reconstruction is made with +1 grating 
order F7 2 (p)R*(p) for p = 1 (zone 1), and p — 2 (zone 2), 
and to the -1 grating order E’Kp)R(p) for p = 2 (zone 3) 
and p = 1 (zone 4). Because the calculations are made 
with two frames hologram, the zero order terms \R(p )\ 2 
cancel, and are not visible on Fig. Efa). 

Note that the angular tilt of the beam splitter BS as 
well as the source point positions F\ and F 2 (see Fig. [l]) 
have been chosen so that the four regions in Fig. [2fa) do 
not overlap, and have sharp edges [2[. From Fig. [2fa), 
we have selected the desired +1 grating orders E 2 R* by 
cropping zone 1 and zone 2 (for p — 1 and 2) and by 
taking the inverse Fourier transform of the cropped zones: 


E 2 R*(x,y,p) = FFT -1 C P [fi(k x ,k y ,p) 


(6) 


where C p is the crop operator for polarisation p. 

Since R(p) is roughly constant with the position x, y , 
the correlations C 2 (£,£') can be then calculated by re¬ 
placing E 2 (x,y,p,t) by [E 2 R*](x, y,p, t) in Eq. [5j The 
statistical average <|C 2 1 2 ) was obtained by first recording 


the sequences of 150 camera frames: /(to)•••^(£ 149 ), at 
times: t n = nAt and At = 100 ms, yielding 75 holo¬ 
gram: i/(£o), i/(t 2 )...i/(ti 48 ), which were used to cal¬ 
culate [E 2 R*](x,y,p, t 2n ) and t 2n /), and then by 

averaging |C 2 (t 2 n, t2n')| 2 f° r all couple of times t 2n ,t 2n / 
with |n — n'\ > 5 and n,n' = 0..74. 

Because of experimental defects, \R(p)\ varies slightly 
with position. This affects the calculation of C 2 (t, t') and 
(|C 2 | 2 ). In order to account for this effect, we measured 
| R(x, y , p) | from our stack of holographic data and we cal¬ 
culated C 2 (£,£') with [E 2 R*](x : y,p, t)/\R(x, y,p)\- This 
correction is about 10% for (|C 2 | 2 ). 


NUMBER OF GEOMETRICAL MODES N' g 

Figure [2] (d) displays the hologram \H \ 2 we got with¬ 
out sample. |i/| 2 exhibits four bright circular zones that 
are smaller in diameter than with the sample (Figj2^a)). 
These circles correspond to the MOl pupil located in the 
plane P’ that appears sharp in the plane P, because MOl 
and M02 form an afocal optical system. There is thus 
no field out of the MOl collection angle, and the number 
of geometrical mode N'g must be calculated with MOl 
numerical aperture NA=0.9. 

We must notice that N' g is a little bit smaller than the 
number of pixels of zones 1 and 2 in Fig. [2] (d). Similarly 
N g is a little bit smaller than the number pixels of zones 
1 and 2 in Fig. [2] (a), but the difference is larger, because 
the brightness within the pupil decreases noticeably near 
the pupil edge. This means that the reconstructed field 
within the pupil is random from one pixel to the next. 

We used this property to calculate N g and TV'. We 
measured the averaged intensities (|/^| 2 ) for each pixels 
of zone 1 and 2, and we used this information to calculate 
by Monte Carlo the residual correlation (|(C 2 | 2 ) that is 
expected for pupils fields random in space and time. The 
number of mode N g and N g we got by this way, with 
and without sample, agree within a few per cent with 
2 tt[NA] 2 L 2 /A 2 with NA=1.4 and 0.9. 


[1] E. Cuche, P. Marquet, and C. Depeursinge. Spatial fil¬ 
tering for zero-order and twin-image elimination in digital 
off-axis holography. Appl. Opt., 39, 4070-4075, 2000. 

[2] N. Verrier, D. Alexandre, G. Tessier, and M. Gross. Holo¬ 
graphic microscopy reconstruction in both object and 
image half-spaces with an undistorted three-dimensional 
grid. Appl. Opt., 54, 4672-4677, 2015. 
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Light can be transmitted through disordered media with high efficiency due to the existence of 
open channels, but quantitative agreement with theory is difficult to obtain because of the large 
number of channels to control. By measuring, by digital holography, the correlations of the transmit¬ 
ted field, we are able to determine the number of channels involved in transmission. With sample, 
the number of channels is divided by about 15 with respect to what is expected without sample. 

This figure is in good agreement with the theoretical prediction. 

PACS numbers: 42.25.Bs Wave propagation, transmission and absorption; 05.60.Cd Classical transport; 

02.10.Yn Matrix theory; 42.40.Ht Hologram recording and readout methods. 


Wave propagation in complex media gives rise to a 
multitude of effects that physicists are trying to under¬ 
stand and analyze both experimentally and theoretically. 
Some of these effects, such as coherent backscattering 
were observed mi , and are well understood. Others, 
like the 3D Anderson localization of an electromagnetic 
wave, are a kind of Holy Grail Q that may exist or not. 
Other, are under study. For example, recent wavefront 
shaping and phase recording technologies make possible 
to focus and image through strongly scattering media 0- 
0 ], or to get enhanced transmission by adjusting the input 
mode structure to fit with the so-called ’’open channels” 
In the best case, the transmission was increased 
by a factor of 4 with respect to that observed without ad¬ 
justment @, while optimal transmission remains small as 
compared to the 100 % expected with open channels [11]. 

Two methods have been used to study this open chan¬ 
nel effect. The first method is to change the mode struc¬ 
ture of the incident field with a spatial light modula¬ 
tor (SLM) to maximize the transmission [10]. Once the 
maximum is reached, it is considered that the incident 
mode fits with an open channel, and that the obtained 
transmission is that of an open channel (ideally 100%). 
The second method is to measure and to diagonalize the 
transmission matrix t 0]. The eigenvectors, whose eigen¬ 
values are maximal, correspond to open channels [8|, 9]. 

The study of open channels in an open geometry (slab) 
by both methods presents great difficulties because of the 
large number N g of geometrical modes (or channels) to 
control or measure. Indeed, the observation of signifi¬ 
cant effects requires a medium, whose thickness / is large 
compared to the light transport mean free path /*. Fur¬ 
thermore, lateral dimensions L should be large compared 
with the thickness l so that the effect can be analyzed 
easily. So we must have L l /*, i.e. for example 
l ~ 10/* and L ~ 100/*. On the other hand, it holds 
N g = 2ttS/X 2 (for 2 polarisations), where S = L 2 is the 
sample area, and A the wavelength. Since /* ^ A in the 
best case, we get N g ~ 10 5 . This number is much higher 
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FIG. 1: Typical distribution of the eigenvalues t n of t t 1 plot¬ 
ted in decreasing order. Eigenvectors are calculated by using 
a typical transmission matrix [t] of size N x N with N = 256 
obtained by Monte Carlo. The parameters of the calculation 
are adjusted to get an average transmission Y^n= i ^n — 0.1. 


than the number of modes controlled or measured in ex¬ 
periments 00. 

Popoff et al. 0 control only 650 modes. They get an 
increase in the transmission by a factor 3, but nothing 
indicates that the optimization process is completed, and 
that the obtained maximum of transmission corresponds 
to an open channel. Yu et al. [9] perform measurements 
on a base of 21078 modes, and the distribution of the 
eigenvalues they get is quite far from that predicted by 
the open channel theory. 

Rather than preparing open channels which is a diffi¬ 
cult task, we propose to study the outgoing field, which 
is directly observable. This field is distributed over a 
number N 2 of modes, which is about the number of open 
channels, that is N g l*/l 0 , this is thus much smaller 
than the number of geometrical modes N g . For random 
incoming fields, one expects thus to observe larger corre¬ 
lations for the outgoing field than for the incoming one. 
In the following it is shown that these correlations can in¬ 
deed be used to define and measure N 2 , that they provide 
a mean to highlight the open channel effect by measuring 
the ratio N g /N 2 1, and to make quantitative compar- 
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ison with theory. 

To be more specific, let us consider the eigenvectors 
U n of t^t where t is the transmission matrix. These 
eigenvectors U n , whose eigenvalues are t 2 (where t n > 0 
is real) constitute a natural orthonormal basis of the in¬ 
coming field Ei with U n .U n / = 5(n,n'). Similarly, the 
eigenvectors V n of t C constitute an orthonormal basis 
of the outgoing field E 2 with V n .V n / = S(n , n'). U n and 
V n are related together by tU n = £ n V n . We calculated, 
by Monte Carlo, a typical transmission matrix t obtained 
for N g = 256 modes in the random matrix theory ideal 
case, and we calculated the eigenvector U n and the eigen¬ 
value 1 2 . Like Choi et al. jl3[, we have plotted on FigJT] 
the eigenvalues t 2 in decreasing order from t 2 = 1 for 
n = 1 down to t 2 ~ 0 for n = 256. The calculation was 
made for l = 9/* to get an average transmission for en¬ 
ergy T = (1 /TV) Yln=i c l° se to 0.1. The open channels 
correspond here to the first eigenvectors, with — 1, 
while the outgoing channels correspond to the channels, 
whose transmission t 2 is noticeable (n = 1 to about 40). 
Ei is distributed over all modes U n with n = 1 to 256, 
while the outgoing field E 2 is mainly distributed over the 
so-called ” outgoing modes” V n with n = 1 to about 40. 

Let us consider a typical input field Ei = Yln =1 c nUn, 
where c n are independent identically distributed Gaus¬ 
sian complex random numbers. The outgoing field is 
E 2 = tEi = En«i c n^nV n . The incoming and out¬ 
going energies are respectively |Ei| 2 = Yln =1 l c n| 2 and 
|E 2 | 2 = J2n =1 ^nl c n| 2 - The Monte Carlo calculation gives 
a transmission efficiency |E 2 | 2 /|Ei| 2 close to T = 0.1. If 
E 2 is re-injected into the sample by using a phase con- 
jugaison mirror, we get the twice transmitted field (i.e. 
the outgoing field that is back transmitted through the 
sample) E 3 = J2n=i c n 2 2 U n , whose energy is |E 3 | 2 = 
Enli^M 2 . We can define the back transmission Tb 
of the outgoing field by = | E 3 1 2 /1E 2 1 2 - We got Tb= 2/3 
as expected m. 

In this letter, we propose to study the ” equivalent num¬ 
bers of modes” Ni and N 2 of the incoming (Ei) and out¬ 
going field (E 2 ). Since Ei is uniformly distributed over 
all the N g geometrical modes U n , one has N\=N g . On 
the other hand, it holds N 2 «C N g since E 2 is distributed 
over a much smaller number of outgoing modes V n . The 
” equivalent number of modes” N 2 corresponds here to 
the ” width” of the curve t 2 (i.e. N 2 —40 for the example 
plotted on Fig|Tj) . 

In order to define N 2 more precisely, we considered two 
uncorrelated random incoming fields Ei = Yln=i c nU n 
and E’i = Ylm=i c n^n having the same energy |Ei| 2 = 
|E’i| 2 , and the corresponding outgoing fields E 2 = 

J2n =1 c ntnV n and E’a = J2n=i c «^V n . One can then 
calculate the normalized residual correlation C\ and C 2 



FIG. 2: Scheme of the experimental setup. 42(1), 42(2): refer¬ 
ence fields of polarisation p—1 and p=2; BS, PBS, beam split¬ 
ter and polarized beam splitter; DL: diffusing liquid; MOl and 
M02: microscope objectives; Q, Q’: camera plane, and cam¬ 
era conjugate plane with MO2; S: sample and sample outgoing 
plane. 


for the incoming and outgoing fields. 
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We calculated the statistical averages (|Ci| 2 ) by Monte 
Carlo and got (|Ci| 2 ) = 1/N g = 1/N\. This result can 
be proved rigorously (see the supplementary material). 
This suggest to define the equivalent number of mode 
N 2 by (IC 2 I 2 ) = 1/AT 2 . We have then calculated typical 
transmission matrices t for 7^=256, 512 and 1024 and 
for ///*=1 to 40. For each t we have calculated (IC 2 I 2 ) 
(and thus 7V 2 ) by Monte Carlo. For l/l* 1, we got the 

simple empirical formula: 


Ci = 

c 2 = 


N 2 = (3/2 )TNg (2) 


The number of outgoing mode N 2 is thus 3/2 x larger 
than number of open channels ( TN g ). 

The key idea of this work is that (IC 2 I 2 ) and 7V 2 can 
be easily measured experimentally, with a time varying 
random incoming field, which is uncorrelated at times 
t and t'. Indeed, one can get the correlation C 2 from 
the measured outgoing field ^(x, y, t,p) (where x,y are 
the transverse coordinates, t the time, and p = 1,2 the 
polarisation) by: 


„ U - E x ,y, p E2(x,y,t,p)EZ(x,y,t',p) ^ 
2(,) " (3) 


The importance of open channels can be directly demon¬ 
strated by measuring 7V 2 = 1 /(|C 2 | 2 ) and by verifying 
that 7V 2 /A^<Cl. Furthermore, a quantitative comparison 
with our empirical model can be made by verifying EqE] 
experimentally. 

Figure [2] shows the scheme of holographic setup. The 
sample S is a ZnO powder slab with thickness l = 
22 /im ± 7 /im deposited on a microscope cover slide. In 
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FIG. 3: (a) Averaged correlation (|Cj 2 ) measured with ZnO 

plotted as function of MO2 position z. Calculation is made 
with M x M pixels with M m 1024 (red diamonds), 362 
(green crosses), 256 (bleue rectangles), 181 (purple crosses), 
128 (dark blue triangles) and 90 (brown stars). Solid grey 
lines is l/^=3.83e-6. (b) Averaged correlation (\C'\ 2 ) mea¬ 
sured without ZnO. Solid grey line is l/iV^=9.286e-6. 


order to maximize the collection of both input and output 
modes, the sample is positioned between two microscope 
objectives: MOl (NA = 0.9 air, x60) in the powder side, 
and M02 (NA = 1.4 oil, x60) in the cover slide side. A 
tank (DL) filled with viscous diffusing liquid is positioned 
in front of MOl to randomize the illumination structure 
in both time and space. The outgoing fields E 2 (x, y , t,p) 
interfere with two orthogonally polarized reference beams 
R(p) whose off axis angles are slightly different. It is then 
possible to extract E 2 (x, y , t,p) from the interference pat¬ 
tern I = i 2 | R(p) + E 2 (p)\ 2 that is recorded by the 
camera. This calculation involves 2 frames holograms 
H — I n — J n+ 1 where I n and / n +i are successive camera 
frames. See supplementary to get details on experiment 
and on data analysis. It should be noted that the exper¬ 
imental tradeoffs (i.e. diffusing liquid randomized illu¬ 
mination instead of SLM and two frame holograms) were 
made in order to minimize the correlation bias within the 
incoming fields E\. In particular, two frames holograms 
minimize bias due to ballistic light that travel through 
the diffusing tank DL. 

To have an extra control parameter, the M02 position 
along z is adjusted by a motorized micrometer device. 
For z ~ 0, the camera’s conjugated plane Q’ coincides 
with the exit plane S of the sample, and plane S is on 
focus on the camera. This configuration maximizes the 
open channel effect. When \z\ ^ 0 increases, the mea¬ 
surement plane Q’ no longer coincides with S and (IC 2 I 2 ) 
decreases. Indeed, if Q’ 7 ^ S, the correlation C 2 is calcu¬ 
lated without being able to define precisely the sample 
selected zone where the outgoing field is measured: points 
outside of the zone contribute to the signal, and points 
inside the zone contribute less since the detection angle is 
lower [16]. Controlling z provides thus an experimental 
mean to decrease and even to suppress the open channel 
effect. This is expected to occur when z becomes much 
larger than the sample transverse size L. In that case, 
(IC 2 I 2 ) must be close to (|Ci| 2 ) = 1 /N g . Data are ac¬ 
quired by moving M02 from z=Q’S=-80 to 250 pm by 


step Az=10 or 20 pm. For each position z of M02, 150 
frames I are recorded and <|C 2 1 2 ) is calculated with 75 
holograms H . 

In order to vary the sample transverse size L, the data 
recorded by the camera (1340 x 1024 pixels) are cropped 
to M x M with L = 77 pm for M = 1024. The average 
correlation (IC 2 I 2 ) is multiplied by M 2 /1024 2 to com¬ 
pensate for area variation. This provides a second ex¬ 
perimental mean to decrease the effect of open channels 
when the transverse size L of the sample becomes smaller 
than its thickness /, since, in that case, the transverse 
geometrical losses becomes large. Conversely, by extrap¬ 
olating L toward infinite, one can get an estimate of the 
maximum enhancement of correlation that is expected 
for a given thickness l of the sample. 

We have verified experimentally that changing z and 
M does not modify the detected energy per unit of area. 
We have plotted on Fig. [3] (a) M 2 (|C2| 2 )/1024 2 as a func¬ 
tion of the M02 position z for M = 1024 (red diamonds), 
362 (green crosses), 256 (bleue rectangles), 181 (purple 
crosses), 128 (dark blue triangles) and 90 (brown stars). 
As expected (|C 2 1 2 ) is maximum for z ~ 0, and decreases 
with \z\. When the sample size M decreases, the geomet¬ 
rical losses increase, and M 2 (\C 2 \ 2 ) /1024? decreases too. 
Moreover, the peak of correlation becomes narrower, and 
the correlation reaches a plateau for small M and large £ 
(e.g. M = 90 and 2 = 250 pm). This plateau is reached 
when geometrical losses are large, that is when z /, L. 

We compared this plateau with the correlation 
<|Ci| 2 ) = 1 /N g that is expected without open channels. 
The usual relation that defines the number of geomet¬ 
rical mode N g should be corrected to account for the 
refractive index n of the medium where the output field 
is collected (factor xn 2 ), and for the detection solid angle 
(factor x[NA] 2 /n 2 ) where NA is the numerical aperture 
of the detection objective. One gets thus: 

N g = 27r[NA] 2 £ 2 /A 2 (4) 

where L is the transverse size of the sample and A 
the wavelength in vacuum. For NA=1.4 and L=77 
pm, we get l/7Vg=3.83e-6 (solid grey line on FigEJa)) 
in good agreement with the plateau of correlation 
curves M 2 (\C 2 \ 2 )/1024 2 . This agreement validates our 
correlation-based theoretical approach and proves that 
we have correctly evaluated the number of geometrical 
modes by EqlH . 

We conducted a control experiment without ZnO sam¬ 
ple and we calculated the correlation (IC^I 2 ) versus 2 for 
M = 1024. As seen on Fig|3] (b), (|C 2 1 2 ) * s mu ch lower 
than with the sample, and does not depend on z position. 
Moreover, (IC^I 2 ) —8.8e-6 is close to the correlation that 
is expected without open channels (|C^| 2 ) = 1/N' g =9.3e- 
6 (see solid grey line on FigEfb)), where N'g is calculated 
by Eq]4]with MOl numerical aperture NA=0.9 (see sup¬ 
plementary) . 
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x=1024 jM 


FIG. 4: Maximum averaged correlation 1024 2 (|C2| 2 )/M 2 ac¬ 
cording to x — 1024/M. Measured point are red diamonds. 
Solid green line is a±e~ x ^ hl + a 2 e~ x ^ b3 + <23, with ai=5.19e-5, 
a 2 =1.59e-5, &i=1.7, b 2 = 10.0 and a 3 =3.67e-6~l/iV 

To calculate the correlation that is expected without 
lateral geometrical losses of you slab sample, we have 
plotted on Fig. 0] (red diamonds), the maximum of 
M 2 (\C 2 \ 2 )/1024 2 obtained during the 2 scan, as a func¬ 
tion of x = 1024/M. We extrapolated our results to 
L = 00 by fitting the data with a sum of 2 decreasing 
exponentials (see solid green line on Figj4]). Extrapola¬ 
tion to L= 00 (i.e. x = 0) yields M 2 (\C 2 \ 2 )/1024 2 =7.lbe- 
5. Half maximum is reached for x=1.73 i.e. L=50 /am 
(about twice the thickness /). The importance of open 
channels that corresponds to N g /N 2 1 is here clearly 

demonstrated, since we got A^/7V 2 =12.3 without extrap¬ 
olation for L — 77 /im, and 18.6 after extrapolation to 
L = 00 . 

The N g /N 2 measured ratio must a priori be cor¬ 
rected by the collection angle of the objectives MOl and 
M02 [14[. M02 losses yield a geometrical loss factor 
[NA] 2 /n 2 =0.85 that affects N g (see Eq0|). Monte Carlo 
calculations Ed show that they affect N 2 by roughly the 
same amount when N g /N 2 1. The ratio N g /N 2 is thus 
not modified by the M02 losses. On the other hand, the 
MOl losses decrease the number of incoming geometri¬ 
cal modes N\ = TV' by a factor 0.9 2 /1.52 2 =0.35. This 
factor was accounted for to get the l/N'g solid grey line 
displayed on Figj3] (b). Calculations show that this de¬ 
crease of N± by 0.35 yields a small decrease of N 2 (~ 0.9) 
as noticed by Kim et al. 0. The corrected extrapolated 
ratio is thus N g /N 2 =1^.7A. 

To perform a quantitative comparison of our results 
with Eq|2] empirical formula, we measured the transmis¬ 
sion T. The ratio of the signal energy (^ |^ 2 1 2 ) detected 
with and without sample is 1/30 d= 20 % [l8l . Tak¬ 
ing into account the M02 losses (x0.85) (l9j, we get 
T = (1/30)/0.85 ~ 1/25. The ratio that is expected 


from Eq|2]is thus N g /N 2 =2/3T=16.6 dk 20 %, in good 
agreement with the corrected extrapolated experimental 
results 16.74. Although this agreement is already ex¬ 
tremely encouraging, the proposed approach can still be 
improved to provide better quantitative results. A more 
precise measurement of the sample thickness Z, and av¬ 
erage transmission T should be made. A setup using a 
high NA oil objective for MOl yielding roughly the same 
number of geometrical modes for input and output (i.e. 
N g ~ N f g ) would also give a better highlight of the open 
channel effect. 

More fundamentally, the shapes of Fig]3] (a) network 
curves, and of FiglU curve give informations that were 
exploited not. Indeed, these curves characterize the de- 
pendance of (|C 2 |) with the M02 defocus z, the sample 
transverse size L and the sample thickness /, i.e. with 
the slab geometry. In order to exploit these informations 
a theory that accounts for the slab geometry should be 
developed. We hope that the present results will en¬ 
courage such development. We acknowledge ANR-2011- 
BS04-017-04 grant for funding. 
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